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Copper NMR has been used as a local probe of the oxygen ordering in Ortho-II YBa2Cu306.5 
crystals grown in BaZrC>3 crucibles. Line assignments have been made to each of the expected 
crystallographically inequivalent sites. The presence of distinct and narrow lines for these sites as 
well as the lack of a line known to be associated with oxygen defects indicates that these crystals 
are highly stoichiometric. Our estimate of the lower limit on the chain length is consistent with that 
derived from X-ray diffraction measurements. In addition, we have found no evidence for static 
magnetic moments, in contrast to some previous results. 



I. INTRODUCTION 

Since the discovery^] of high temperature supercon- 
ductors (HTSC), considerable effort has been devoted 
to understanding the electronic and magnetic proper- 
ties of these materials at a microscopic level. Although 
there are a large number of perovskite HTSC materi- 
als, the YBa2Cu3C>6+a; (YBCO6+2;) family has a special 
place. It was the first superconductor discovered|2j with 
a transition temperature T c above liquid nitrogen tem- 
perature with a relatively simple synthesis route. It is 
available in almost the entire region of the phase dia- 
gram, easily tuned by varying the oxygen content in a 
quasi-continuous way. Thus it is appealing for probing 
the intrinsic properties of HTSCs. In the underdoped 
regime, a variety of crossover phenomena are observed at 
temperatures above T c in which various forms of spectral 
weight at low energies are suppressed. These phenomena 
are associated with the opening of a pseudogapQ. To 
help elucidate the nature of this pseudogap, it is desir- 
able to have highly stoichiometric crystals with a high 
degree of crystalline order. In particular, such samples 
are needed if one wants to determine whether or not the 
crossover behaviors are transistions smeared by disorder 
and/or impurities^- 

YBC06+X can be prepared in the whole range of oxy- 
gen content 6 + x = 6.0 to slightly larger than 7.0. The 
oxygen content affects the crystallographic structure of 
the material 5]. At £=1.0, the crystal structure con- 
sists of a bilayer of two-dimensional (2D) Cu02 planes 
and a single layer of one-dimensional (ID) Cu-O-Cu 
chains, hereafter denoted the CuO chain layer. Nuclear 
magnetic resonance (NMR) spectroscopy offers a unique 
probe to investigate the microscopic properties of the 
Cu02 planes and CuO chain layers separately. In par- 
ticular, the response of inequivalent sites can in principle 
be resolved in NMR spectra and the dynamical prop- 
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erties of these sites studied by nuclear relaxation. The 
chains are important not only because of their role as 
charge reservoirs and their proximity to the supercon- 
ducting Cu02 planes, but also because of their low di- 
mensionality. It is well known that charge/spin density 
wave(CDW/SDW) ground states are prevalent in quasi- 
one dimensional systems 0. Therefore, the chains in this 
system present an important opportunity for clarifying 
the question of whether such transitions exist in HTSC. 
In fact this has been the subject of study in the stoichio- 
metric compound YBa2Cu40s containing two CuO chain 
layers, where evidence for the presence of incommensu- 
rate CDW fluctuations, near the wave vector 2kp occur- 
ring below 180 K has been reported^. There are a few 
reports on charge modulations in YBC07 as well S]. Per- 
haps the most striking evidence for charge oscillations in 
the chains of YBC07 is provided by STM observations 0- 

The structural phase diagram of YBCO6+3; has 
been studied in great detail; early neutron diffraction 
reveal^^ a tetragonal to orthorhombic phase transition 
as the sample is cooled from the growth temperature, 
in which oxygens form ID CuO chains and thus break 
the tetragonal symmetry. Typically the chain ordering 
occurs in domains and the samples need to be mechani- 
cally detwinned in order to produce crystals with a single 
orientation of chains. 

Depending on the preparation conditions and oxygen 
content, oxygen atoms arrange themselves in different 
types of ordering in the CuO chain layers. Efforts to 
produce crystals with domains of long chains requires 
an understanding of oxygen ordering in the chain layer. 
In addition to the orthorhombic structure (Ortho-I, full 
chains) at x—1.0 (Fig. ^ (ii)) and tetragonal phase (T, 
empty chains) at x=0.0 (Fig. H](iii)), at least four more 
modified orthorhombic structures have been observed|llj 
such as Ortho-II (Fig. □ (i)) for 6.35 < 6+x < 6.62 
with alternating full and empty chains. There have been 
theoretical studies based on lattice gas models to ex- 
plain the structural phase diagram of YBC06+x[l3. 
These models can be shown to be equivalent to the Ising 
model with asymmetric next nearest neighbor interac- 



tions (ASYNNNI model) which use three pair potentials 
between oxygen atoms in CuO chain layers: V\ for near- 
est oxygen neighbor pairs, V2 for next nearest pairs con- 
nected through a Cu ion and V3 for those without Cu 
in between, with the condition V% < < V3 < V\. It is 
the strongly repulsive interaction Vi of two inter-chain 
nearest neighbor oxygen ions which prevents the occupa- 
tion of nearest neighbor oxygen sites and promotes the 
formation of chains. Long chain segments are induced by 
the attractive potential V2 between two intra-chain next 
nearest neighbor oxygen ions connected through a Cu 
ion. Alternating sequences of full and empty chains fol- 
low from the weak repulsive V3 between two inter-chain 
next nearest neighbor oxygen ions. 

In YBC06+X, T c changes from the maximum 92.5 
K at optimal doping to zero for fully oxygen depleted 
samples This change takes place in a nonlinear man- 
ner with a pronounced plateau at 60 K for 6.45 < 6+x < 
6.65 and a broad maximum at 92.5 K for 6.85 < 6+x < 
7.0. These features are associated with the Ortho-II and 
Ortho-I ordered phases, respectively. It is suggested 0] 
that this correlation has to do with charge transfer rather 
than with any specific T c of the ordered phases. The 
importance of oxygen ordering in the CuO chain layers 
in determining the charge transfer was first noted in a 
study of the time dependence of T c du ring room tem- 
perature annealing of quenched samples With the 
oxygen content (average occupancy of different oxygen 
sites) remaining constant, the degree of oxygen ordering 
increased with time. This resulted in an increase in the 
orthorhombicity and in T c . It was shown that the re- 
arrangement of oxygen ions in the chain layers to form 
longer chain segments was responsible for this increase. 
Another experimental indication of the presence of a min- 
imum chain length for the charge transfer to take place is 
the Tjv vs. oxygen content for small x: Xjv is nearly con- 
stant for 6+x <6.2 0| . The relationship between oxygen 
ordering and hole concentration has been addressed|l7j 
theoretically by band structure calculations. Chain frag- 
ments were found to be unable to transfer charge to the 
planes until they exceeded a critical length equal to three 
(containing three Cu and two O atoms). An alternative 
model combining the ASYNNNI and the charge trans- 
fer models has also been|l8| successful in explaining the 
variation of T c vs. x. Here it is assumed that only the 
ordered phase participates in superconductivity. Thus 
there is a minimum domain size below which no charge 
transfer can take place. 

NMR experiments can be used to measure the electric 
field gradient (EFG) in materials. Since the EFG at a nu- 
cleus depends on the symmetry of that site, the valence of 
the ion itself and the charge of its environment, valuable 
information about the electronic structure of the system 
under study can be obtained by such measurements; i.e. 
charge transfer. However, the interpretation can be com- 
plicated if the samples are not single phase or suffer from 
impurities and/or disorder. There are a number of NMR 
reports on underdoped YBC06+X with different oxygen 
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FIG. 1: The crystal structure of YBCO6+2 for (i) x=0.5 
(Ortho-II) (ii) x=1.0 (Ortho-I) and (iii) x=0.0 (T). 



contents 0, |2jJ, El H2, Hj| where either the sample is 
not stoichiometric or suffers from oxygen inhomogeneity, 
powder averaging or twinning. Liitgemeier's group |22j 
was the first to study the effects of oxygen ordering on 
the NMR properties of YBCO6+2;. However, they used 
aligned powdered samples where the NMR resonances 
were broadened by non-stoichiometry and oxygen dis- 
order. For Ortho-II ordering only a partially resolved 
splitting of the resonance line from the Cu02 plane was 
observed. 

With recent enhancements in crystal growth very 
pure and highly ordered detwinncd Ortho-II crystals 
with lar ge oxygen-ordered domain size have become 
available |24|. To probe the intrinsic properties of the 
Ortho-II phase of YBCO we have performed NMR ex- 
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periments on highly ordered detwinned Ortho-II single 
crystals of YBa 2 Cu 3 6 . 5 (YBC06.5). We report NMR 
spectra and line assignments for these Ortho-II crystals 
together with NMR evidence pertaining to the quality of 
the crystals, specifically regarding the chain lengths and 
the lack of oxygen disorder. We also report EFG param- 
eters for different Cu sites in Ortho-II structure along 
with an analysis regarding lattice and valence charge con- 
tributions. We will discuss the relation between oxygen 
ordering and charge transfer in the Ortho-II structure. 
It will be shown that these crystals are an extremely im- 
portant disorder-free prototype of the underdoped state 
of HTSC. 



II. EXPERIMENT 

Single crystals of YBC06.5 were grown by a flux 
method using BaZrOs crucibles p4j|. The crystals were 
first detwinned and then the oxygen content set to 6.50 
during the final annealing process. The NMR properties 
of two such crystals used in our experiments proved to be 
identical to within experimental error. Hence we will not 
distinguish between them. The crystals were 2 mm x 2.5 
mm wide and ~ 0.1 mm thick platelets with the c-axis 
perpendicular to the platelet and a and b axes parallel 
to the edges of the platelet. The superconducting transi- 
tions of the crystals measured by a SQUID magnetome- 
ter give T c = 62 K with a transition width AT C = 0.6 K. 
X-ray diffraction measurements performed on these crys- 
tals show a three-dimensional Ortho-II oxygen ordering 
with long correlation lengths 148, 430 and 58 A along the 
a, b and c crystallographic axes, respectively^. The 
long correlation length along the b (chain) axis indicates 
the presence of much longer chains than have been avail- 
able previously 22J. Both the large oxygen-ordered do- 
main size and a small AT C are indicative of high degree 
of oxygen homogeneity in these crystals. Because twin 
boundaries promote oxygen clustering, the removal of 
twin boundaries was found to be essential for obtaining a 
highly ordered Ortho-II phase[24|. Single crystal X-ray 
diffraction measurements on our crystals showed at most 
10% of the sample was twinned. 

Cu NMR data were taken between 30 and 250 K and 
magnetic fields up to 9 T with a home built pulsed spec- 
trometer with quadrature detection. The crystal was 
placed in either a Pt or Au rectangular coil tunable from 
30 to 100 MHz. A rectangular coil was chosen to conform 
to the shape of the crystal and so maximize the filling fac- 
tor. The onset of the superconducting transition of the 
sample was observed in situ at 62 K by a sharp change in 
the inductance of the tuned coil at zero applied magnetic 
field. A powder Al reference sample was located in the 
coil adjacent to the sample to calibrate the value of the 
magnetic field at the sample position. 

The data were obtained with a standard two-pulse 
spin-echo sequence, with 4 /is pulse length and 30 to 
40 /is pulse separation. Copper nuclei have two isotopes 



63 Cu and 65 Cu with spin |. Each isotope gives one 
central transition (| <-* — h) an d two satellites (±| <-> 
±|). To obtain the complete Cu NMR spectrum of the 
sample we used either field or frequency sweep methods. 
The spectra were then constructed by superimposing the 
Fourier transform spectra of the spin echo measured at a 
certain frequency or field interval. There was no attempt 
to correct the signal intensity for Ti or frequency depen- 
dencies in frequency sweep measurements, given the lat- 
ter's narrow range. Although the samples are high qual- 
ity, because of their small volume a significant amount 
of signal averaging was required to obtain a satisfactory 
signal to noise (S/N) ratio. 

III. RESULTS AND DISCUSSION 
A. Line Assignment 

Typical 63 Cu NMR spectra of central transitions ob- 
tained at 60 K at 75.75 MHz for the applied field parallel 
to the a-axis (Ho | |a) are shown in Fig. [21 Similar spectra 
were also observed for Ho||b containing four resonance 
lines but at different positions. We denote the observed 
lines in Fig.[3by A, B, C and D. The spectrum for A-D 
resonances has also been obtained for 65 Cu with inten- 
sities consistent with the natural abundances of the two 
isotopes. These four Cu resonances are associated with 
the four inequivalent Cu sites in the unit cell as expected 
from the Ortho-II crystal structure (Fig.^i)). Two lines 
A and C appear nearly with the same width, integrated 
area and relaxation rates. Line B has a very long spin- 
lattice relaxation time (nearly three order of magnitude 
longer than A and C) and is relatively narrow compared 
to the rest of the lines. Line D is broader compared to 
line A and C and has an order of magnitude longer re- 
laxation time. It should be noted that these resonance 
lines are well resolved and that there appears to be no 
background signal under these lines. Also note that a rel- 
atively long repetition time of 250 msec must be used in 
order not to saturate line B. Most of the spectra are col- 
lected with a shorter repetition time of 25 msec in order 
to saturate line B which otherwise obscures line A. 

Line assignments will be made by comparing the local 
parameters of the four Cu sites with those of YBC06 and 
YBC07. To motivate this comparison consider the struc- 
tural environmemt of each Cu site. Ortho-II YBCO con- 
sists of alternating full and empty chains (see Fig. ^i)) 
and hence there are two inequivalent Cu sites in CuO 
chain layers, one in full Cu(lF) and one in empty Cu(lE) 
chains. Correspondingly there are two inequivalent Cu 
sites in the Cu02 planes, adjacent to full and empty 
chains Cu(2F) and Cu(2E), respectively. Comparing 
the crystal structure of Ortho-II with that of YBC07 
(Fig. n^ii)), one concludes that the local structural en- 
vironments of Cu(2F) and Cu(lF) are similar to that of 
Cu(2) in the plane and four-fold coordinated Cu in the 
chain layers Cu(l)4 in YBC07, respectively. Similarly, 
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FIG. 2: M Cu NMR spectra for central transitions at 60 K and 
75.75 MHz for Hoj|a taken at different repetition time, T rep 
(top panel: T rep =250 msec, middle panel: T rep =25 msec, and 
lower panel: the difference). The horizontal arrow indicates 
the possible range for Cu(l)3. See text for details. The two 
narrow and small signals in the far right side of the spectra 
are from Cu metal in the sample probe. 



the Cu(2E) and Cu(lE) sites have the same local struc- 
tural environment as Cu(2) in the plane and two-fold 
coordinated Cu in the chain layers Cu(l)2 in YBC06. 
NMR parameters for these lines will be compared with 
the known [H |H ^ NMR parameters (Knight shifts, 
nuclear quadrupole frequency and asymmetry parame- 
ters) of the corresponding sites in YBC07 and YBC06 
to determine the line assignments. 

In order to obtain the NMR parameters for the ob- 
served lines, we have measured the complete 63 / 65 Cu 
spectrum including the satellite and central transitions 
at different frequencies and fields along the three crystal- 
lographic axes. A typical complete spectrum is shown in 
Fig.|21for Ho | |c at 70 K. Note that for this direction of the 
applied field, the central and lower satellite transitions of 
lines A and C overlap. The spectra were analyzed us- 
ing second order perturbation theory for the quadrupo- 
lar shift to the Zeeman interaction (see Appendix) for 
central transitions and by exact diagonalization for Ho 
parallel to crystallographic axes for the satellites since 
second order perturbation theory is insufficient for these 
resonances. 

If one writes the resonance frequency v (see Appendix) 
of the central transition to the second order as 
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FIG. 3: The complete Cu NMR spectrum at 70 K and 64.85 
MHz for Ho||c. For this direction, the central transitions and 
lower satellites of A and C lines overlap. 



it is easily seen that K e f / is a linear function of r- Ho y± 
with the intercept equal to K and the slope proportional 
to the quadrupolar coupling constant for that direction. 
For example, in Fig. 0] we have shown this data for the 
central transition of the line D obtained at 70 K and for 
H ||b. We have used the position of the center of grav- 
ity of this line to obtain K e f f since this line is broad at 
low temperatures. The Knight shifts and EFG parame- 
ters are obtained for all sites with a simultaneous fit of 
the central transition and satellite data in the three field 
orientations |28j. 



0.03 



0.02 



0.01 



0.00 




intercept = 0.00234 ± .00025 
slope = 35 ±1 (MHz) 2 



0.0000 



0.0002 



0.0004 



0.0006 



(7 n H ) 2 

transition of the D line at 70 K, Ho lib. 



FIG. 4: K e ff as a function of — „ vi for the 63 Cu central 



The NMR parameters for the four observed lines in 
the normal state are reported in Table [I] For comparison 
we have listed the NMR parameters of Cu(2) and Cu(l)4 
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sites in YBC07 and Cu(2) and Cu(l) 2 sites in YBC06 in 
Table ITU 2a l2a |27| . We can readily assign the four reso- 
nances observed in Ortho-II to the four different crystal- 
lographic sites by comparing their NMR parameters to 
those of Cu sites in YBC07 and YBC06. This compar- 
ison shows that the two lines A and C, with very similar 
features, come from the planes; A from the planar Cu 
sites adjacent to the full chains Cu(2F) and C from the 
planar Cu sites adjacent to the empty chains Cu(2E). 
The EFG principal axis for both sites is along the c- 
axis which is the same for the Cu(2) sites in YBC07 and 
YBC06. The asymmetry parameter for these sites is very 
small (nearly axially symmetric) as with the Cu(2) sites 
in YBC07 and YBC06. The slightly larger asymmetry 
parameter observed in Ortho-II might be explained by 
crystal structure arguments since Ortho-II structure has 
a lower symmetry than the Ortho-I or tetragonal phases. 
The Cu(2F) and Cu(2E) resonances appear at different 
positions because of their different quadrupolar coupling 
constants. Since the quadrupolar coupling depends on 
the ion's surrounding charge distribution as well as the 
local charge at the position of the ion, a different cou- 
pling is expected for the two planar Cu sites as observed. 
We will discuss the contribution of the lattice and valence 
electrons to this coupling in more detail below. The vq 
for Cu(2E) is higher than the observed one in YBC06, 
consistent with the trend of vq vs. doping 29] . The same 
intensity of the Cu(2F) and Cu(2E) lines 30] within our 
S/N is an indication that there are equal numbers of these 
sites in the crystal. Hence there is a high degree of order- 
ing in the chain layers with a periodic sequence of full and 
empty chains. The mere fact that we have observed only 
two distinct Cu(2) signals for both satellite and central 
transitions (Fig. is suggestive of a complete ordering 
of the full and empty chains. In a disordered structure 
the signal would merge to one broad resonance with some 
splitting depending on the degree of ordering. Such par- 
tially resolved splittings for the planar Cu resonances is 
typically in previous generations of YBC06+X samples. 

The resonance line B appears in the same frequency 
range as Cu(2F) but can be distinguished by its much 
longer relaxation time. In fact, among the four observed 
resonance lines, B has the smallest width, Knight shift 
and the longest relaxation time (of order seconds) com- 
parable with that of Cu(l)2 in the insulating tetragonal 
phase of YBC06[3l|. A comparison of the NMR param- 
eters of this line with Cu(l) 2 present in YBC06 with 
empty chains leads us to assign this line to Cu(lE) in 
the Ortho-II structure. Similar to Cu(l)2, the Cu(lE) 
resonance has axial symmetry. The long spin-lattice re- 
laxation time is interpreted as a result of a low conduction 
electron density at the resonant nucleus. This is consis- 
tent with the Cu ions being in the 3d 10 electronic con- 
figuration. A small Knight shift is also consistent with a 
low density of conduction electrons at this site. In addi- 
tion, the fact that the resonance signal from this site is 
very narrow is in agreement with this conclusion because 
insulating sites tend to be much more uniform in nature 



than conducting ones as they lack an inhomogeneous K s 
component. The Cu(lE) linewidth is 40 kHz at 60 K with 
Ho||a = 6.5 T for different directions of applied field and 
becomes even narrower at higher temperatures. 

Line D is broader than other resonances, at least at 
low temperatures. We have found a large value for its 
asymmetry parameter (see Table . The highly asym- 
metric behavior of the EFG is reminiscent of Cu(l)4 in 
YBC07. This fact together with its similar vq has led 
us to assign this line to the Cu(lF) site. An examination 
of Cu(lF) site symmetry explains its highly asymmetric 
nature. In the Ortho-II structure, the Cu(lF) sites are 
surrounded by four oxygen ions, two chain oxygens and 
two apical oxygens. Since these oxygen neighbors are not 
all equivalent, a non-symmetric EFG is expected. To un- 
derstand the nature of the broadening of this line at low 
temperatures, we have studied the line shape and relax- 
ation rates for several directions of the applied field. We 
attribute this broadening to the presence of Friedel-like 
oscillations nucleated at the chain ends. The details of 
this study will be presented elsewhere 32]. 

We have confirmed that the ratio of quadrupolar shifts 
of the two isotopes measured at the same site differs by 
their quadrupole moment ( 63 Q/ 65 Q = 1.079|2g|) for all 
sites. This indicates that no internal magnetic field, due 
to static magnetic moments, is present at these Cu sites 
since such a field would shift or broaden the lines. This is 
in disagreement with some neutron-scattering measure- 
ments where a moment of ~ 0.05 /is is observed|33 in 
YBC06.5 and of - 0.02 /ig in YBC06.6[3i| for cop- 
per sites in the Cu0_9 planes, and with /zSR results on 
YBC06.67 showing (351] the existence of static magnetism 
in these planes. A magnetic moment of the order of ~ 
0.02 [is would have shifted or broadened the Cu(2F/E) 
lines by about 3 MHz which we would have easily ob- 
served. More recent neutron scattering experiments on 
Ortho-II found 36] no evidence for the presence of static 
magnetic moments. Also, in a further work, it was 
shown 37] that the /iSR results might arise from charge 
inhomogeneity rather than magnetic order. 

To further elucidate this apparent discrepency several 
factors should be noted. As reported by Sidis et aLp^ 
their static moments fluctuate on a nanosecond time scale 
which would average out these moments on the NMR 
time scale. However this still does not settle the discrep- 
ancy with the work of Stock et a/.|36|. In this experi- 
ment the Ortho-II sample is much more stoichiometric 
than that of Sidis et al.au the YBC06.6 sample of Mook 
et ai. Thus, as suggested by Stock et ai, oxygen dis- 
order in the chains may be responsible for establishing 
magnetic order there. With respect to the /LtSR results, 
Sonier et al. 37] have suggested that for samples with x > 
6.5 local charge inhomogeneity, caused by CDW in the 
chains as shown by Grevin et ai.@] could lead to mag- 
netic moments via the stripe mechanism. However for the 
Ortho-II concentration of x = 6.5 no /iSR signature of a 
CDW was detected. From our work on the Cu(lF) line 
mentioned above|33 we have established that a Friedel- 
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TABLE I: Experimental NMR parameters for the four observed lines at 70 K. Knight shifts are reported in units of percent 
and vq in MHz. EFG-PA represents the direction of V zz . Site assignment is based on the comparison of these parameters with 
those of Cu sites in YBC07 and YBC06. For details see the text. 



resonance 


site 


K a 


Kf, 


K c 


V Q 


V 


EFG-PA 


A 


Cu(2F) 


0.25±0.01 


0.25±0.01 


1.18±0.01 


30.55±0.03 


0.018±0.002 


c 


B 


Cu(lE) 


-0.4±0.1 


0.07±0.05 


-0.1± 0.1 


32.1 ±0.2 


0.00±0.01 


c 


C 


Cu(2E) 


0.26±0.01 


0.24±0.01 


1.14±0.01 


27.43±0.03 


0.022±0.002 


c 


D 


Cu(lF) 


0.9±0.1 


0.23±0.03 


0.23±0.05 


19.6 ±0.1 


0.8±0.1 


a 



TABLE II: NMR parameters for Cu(2) and Cu(l) 4 sites in YBC07 and Cu(2) and Cu(l) 2 sites in YBC06 materials. Knight 
shifts are in units of percent, vq in MHz. EFG-PA represents the direction of V zz . 

site 1^ Kb Kc vq 77 EFG-PA 

YBC07 Cu(2) 0.59±0.04 0.59±0.04 1.267±0.001 31.54±0.05 0.01±0.01 c 

Cu(l) 4 1.38±0.07 0.55±0.07 0.60±0.04 19.35±0.08 0.95 ±0.02 a 

YBCO6 6 Cu(2) 0.47±0.02 1.335±0.005 23.80 ±0.2 c 

Cu(l) 2 0.18 ±0.01 -0.117±0.003 29.54 ±0.02 c 

"The Knight shifts are from Ref. 25], v Q and 77 are taken from Rcf. 26] all at 100 K. 'Taken from Ref. 27] at 505 K. 



like spin-density oscillation exists in the chains of suffi- 
cient magnitude to produce the field, of order 1 gauss, 
necessary to account for the Ortho-II /xSR results. Thus 
the CDW/stripe mechanism need not be invoked in the 
case of stoichiometric Ortho-II YBCO. 



B. Chain Length 

We now estimate the chain lengths from our NMR 
data. In YBC06±x compounds with intermediate oxy- 
gen content, a resonance line with vnqr — 23 MHz is 
observedj^ and is associated with three-fold coordi- 
nated Cu(l)3 sites corresponding to the terminal Cu ion 
of the full chains with a relative abundance that dimin- 
ishes with increasing length of the chain segments. A 
similar value of vq for Cu(l)3 and Cu(lF) sites together 
with the fact that Cu(lF) is very broad at low tempera- 
tures, makes it difficult to obtain a precise value for the 
concentration of these sites in our sample. As discussed 
above we expect the concentration of Cu(l)3 to be low. 
In fact, in the complete NMR spectrum for Ortho-II mea- 
sured between 1 to 9 T, to within our S/N we did not 
observe any signal that could not be assigned to the four 
inequivalent sites in the structure. Hence, we can only 
estimate an upper limit on the number of Cu(l)3 sites 
from our S/N and ultimately a lower limit on the chain 
length in our sample. 

To put a quantitative limit on the concentration of 
Cu(l)3 sites, we consider two sets of spectra. The first 
one was taken at 250 K where Cu(lF) is considerably 
narrower than the one at 70 K. The Friedel-like oscil- 
lations are thermally smeared out at this high tempera- 
ture. In Fig. [3]we have shown this spectra for Ho||c. We 
have used the established NMR parameters of the Cu(l)3 



sites [HHH {v NQR = 23.77 MHz, rj = 0.3 and V zz \\b) to 
calculate its possible position. The horizontal arrow in 
Fig. [5] covers the range consistent with parameter un- 
certainties. From the relative intensities of Cu(lF), I4, 
and Cu(l)3, J3, one is able in principle to find the mean 
number of oxygen atoms in the full chains (chain length) 
according to the relation 22] : n— 1+214/13. Since we have 
not observed any additional peaks due to Cu(l)3, we can 
only put a limit on the maximum intensity of ^3 from our 
S/N. To obtain this limit, we have assumed the Cu(l)3 
line has the same NQR width as observed by Yasuoka et 
al. [2(| leading to an NMR linewidth of 66 kHz. If we 
assume we can detect such a line with S /N of two then we 
obtain a lower limit of 32 lattice constants for the chain 
length. The inset to Fig. \5\ illustrates the Cu(l)3 signal 
expected for n=30 and n=60. 

The second set of spectra was taken at 60 K with Ho | |a. 
In this case we use the fact that the relaxation rate of 
the Cu(l)3 site is relatively slow.39f to extract only the 
slowly relaxing components. As illustrated in Fig. [21 two 
spectra taken at a repetition time of 25 msec, and 250 
msec, have been subtracted to reveal the Cu(lE) line as 
expected and no other line. As above we put an upper 
limit on the chain length given our S/N. The calculated 
range for Cu(l)3 is displayed in Fig. |21 Our estimate of 
the chain length from this analysis is 50 lattice constants. 

These estimates are consistent with the XRD result 
of 110 lattice constants for the domain size along the 
chain direction. Using the ASYNNNI model we have 
performed Monte Carlo calculations to determine what 
the domain interface looks like for the two cases of (a) the 
case of a twin boundary and (b) parallel chains that are 
out of phase, i.e. FEFE adjoining EFEF (see Fig. EJ. In 
both cases a linear boundary is preferred as the ground 
state with Cu(l)3 atoms at the domain boundary. Thus 
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we expect to have three-fold coordinated copper sites at 
least along the domain boundaries normal to b and more 
along the twin boundaries to the extent that they exist in 
our crystals. Hence the relative concentration of Cu(l)3 
sites should be I3/I4 ~ 2% from the XRD results whereas 
our detection limit for I3/I4 ~ 4%. From the two values 
for n obtained above, the consistent estimate for the lower 
limit of the chain lengths is 50 lattice constants. 
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FIG. 5: 63 Cu NMR spectrum at 250 K and H (= 7.96T)||c. 
The arrow shows the estimated range of the Cu(l)3 resonance. 
Inset illustrates presence of Cu(l)3 line synthetically added to 
the data for two chain lengths. This line is barely visible for 
n=30 and undetectable for n=60. 



-o o 

©- 

O 

« a 

o o Cu(l), ° & 

G- 

o o 



provements in sample preparation techniques have re- 
sulted in the observation of narrower lines ^(j empha- 
sizing the importance of sample quality in determining 
the properties of HTSC materials. In general, the ob- 
served linewidths can be attributed to the distributions 
in either the EFG parameters (arising from the structure) 
and/or the Knight shift. Our analysis of the linewidths 
of the central transitions vs. applied magnetic field as 
well as a comparison between the linewidths of central 
and satellite transitions indicates that the line broaden- 
ing for Cu(2F/E) and Cu(lE) sites is mainly caused by 
a distribution of EFG parameters while for Cu(lF) it is 
mainly due to a distribution of the Knight shift. The 
linewidths of all resonance lines decreases with increas- 
ing temperature. The planar Cu sites for Ho||c have a 
linewidth of 40 kHz at 250 K. To our knowledge these 
NMR linewidths are the smallest observed for YBC06.5 
samples |40j, another indication of high sample quality. 

Oxygen disorder is a feature of non-stoichiometric 
samples. Erb et a/.^(J have compared the NMR spectra 
of YBC07 and YBC06.9. The former has very narrow 
NMR lines. For example the linewidth of Cu(l)4 is 22 
kHz with Ho = 9 T ||c at room temperature. When the 
oxygen concentration is reduced by 0.1 the linewidth in- 
creases to 183 kHz, as a result of oxygen disorder. Also 
noted was the presence of Cu(l)2 with a linewidth of 
13 kHz. They argue that this narrow linewidth is a re- 
sult of oxygen clustering which provides a more uniform 
environment for the Cu(l)2 sites than would be the case 
for randomly distributed oxygen defects. They also argue 
that the absence of a Cu(l)3 signal is consistent with oxy- 
gen clustering as this three-fold coordinated site forms 
the boundary between full and empty chains. As our 
linewidths are similar to that of YBC07 we can infer that 
there are few oxygen defects in our Ortho-II YBC06.5 
crystals. We find that our Cu(lE) linewidth for fields 
along a and b is 30 kHz at 200 K consistent with few 
oxygen defects. Our narrow Cu(lF) line at these high 
temperatures is also consistent with few oxygen defects. 
This is relevant to the Cu(l)3 concentration. Although 
we expect no substantial number of oxygen defects in 
stoichiometric YBC06.5 the absence of a Cu(l)3 line in- 
dicates that our crystals are indeed substantially stoi- 
chiometric. 



(a) (b) 

FIG. 6: Domain interfaces calculated by Monte Carlo cal- 
culations: (a) twin boundary and (b) parallel chains out of 
phase. Line segments represent Cu-O-Cu chains. Chains are 
terminated with Cu(l)3 sites denoted with circles. 

In NMR studies on HTSCs, large linewidths have been 
an obstacle to investigating intrinsic properties. The 
problem arises from the fact it is not readily obvious 
whether these large linewidths arise from imperfections 
of the structure such as oxygen disorder or are an in- 
trinsic property of these materials. However, recent im- 



From the NMR spectrum obtained for Ho||a and b, 
we are also able to investigate the degree of twinning in 
our crystal. Fig. [7\ shows the 63 Cu NMR spectrum for 
Cu(2F) and Cu(2E) at 60 K, 75.75 MHz. A remnant of 
Cu(lE) is observed on the high side of Cu(2F) as our 
pulse sequence to saturate the Cu(lE) signal is not 100% 
efficient. It is evident from the lack of asymmetry of the 
two lines that there is no obvious indication of the pres- 
ence of twinning in our crystal. Quantitatively we find 
by adding an admixture of the twinned spectrum that to 
within our S/N, there could be a maximum of ten per- 
cent twinning present in our sample. This is consistent 
with our XRD measurements discussed above. 
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FIG. 7: 63 Cu NMR spectrum at 60 K and 75.75 MHz for 
Ho||a and b. The spectra are normalized to the integrated 
intensity of Cu(2E) line. Note that the feature on the right 
side of Cu(2F) for Ho||a is due to the imperfect subtraction 
of Cu(lE) signal. This of course does not affect the Cu(2E) 
lineshape. 

C. Charge Transfer and Oxygen Ordering 

Oxygen ordering in the chain layer may have an effect 
on charge transfer to the planes. We now investigate this 
possibility through the planar copper EFG values. Con- 
sidering the intinerant holes in the Cu02 planes, vq for 
the two sites Cu(2E) and Cu(2F) in these planes might 
be expected to be the same. However, we have observed 
a splitting of the Cu(2) resonance. The coupling of the 
nuclei's quadrupole moment with its environment deter- 
mines vq which is directly proportional to the electric 
field gradient at the nucleus. The local EFGs reflect the 
structural and crystallographic aspects of the compound 
and provide information about ionicity and bonding. The 
EFG at the nuclear site is determined by two contribu- 
tions: one arising from the charges on the neighboring 
ions (lattice contribution) V^"' and the second arising 
from the incomplete electronic shells (valence contribu- 
tion) V zz l of the ion under consideration. The total EFG 
may be written as|4l) 

V zz = (1 - loo )V l z f + (1 - R)V™ 1 (3.2) 

where and R are the anitishielding factors which ac- 
count for the antishielding effect of the core electrons of 
the central ion on the lattice and valence contributions, 
respectively. From this equation, if one calculates the 
lattice contribution, one can obtain the valence contri- 
bution from the data. To what extent does the Cu(2) 
splitting stem from chain ordering in the Ortho-II struc- 
ture (lattice contribution) or from differing charge trans- 
fer (valence contribution)? To address this question, we 
have used the point charge model (PCM) to calculate the 



TABLE III: Calculated lattice EFG values, V lat in units of 
MHz, based on the point charge model. 



site 


\rlat 
* aa 


\rlat 
V(,6 


\rlat 
* cc 


rf at 


EFG-PA 


(l-R)V^ 


Cu(2F) 


-1.140 


-1.094 


2.232 


0.020 


c 


-71.62 


Cu(2E) 


-1.204 


-1.078 


2.282 


0.055 


c 


-69.43 


Cu(lF) 


4.182 


-1.791 


-2.391 


0.143 


a 


-96.55 


Cu(lE) 


2.338 


2.476 


-4.814 


0.029 


c 


-2.25 



lattice contributions to the observed EFG values for the 
four Cu sites in Ortho-II structure. Although the PCM 
is simplistic we are interested only in relative differences 
and not the absolute values of EFG parameters. 

Applying the formal valence description to Ortho-II 
YBCO results in Y +3 BaJ 2 Cu^ 2 0^g ionic configuration. 
However this neglects the fact that Cu(lE) is in the 3d 10 
state with charge +1. The resulting extra hole is taken 
to be distributed over the Cu02 planes. We have taken 
the planar Cu sites to have a valence of +2.083, the pla- 
nar O sites a valence of -1.917, the Cu sites in the full 
chain a valence of +2, the Cu in the empty chain a va- 
lence of +1, the rest of the oxygens a valence of -2 and 
Ba and Y in +2 and +3, respectively. We have used 
the atomic positions obtained 42J by single crystal X-ray 
diffraction measurements at 80 K. It should be noticed 
that due to the alternating nature of the full and empty 
chains, some atomic displacements are observed |42j from 
the mean position of atoms, most notably Ba. This is 
attributed to the lattice tendency to minimize its electro- 
static energy [43^. We have summed to convergence over 
the contributions within a crystal containing ±50, ±100 
and ±30 unit cells along the a, b, and c-axes, respec- 
tively. The resulting values for the principal components 
of EFG for all sites are given in Table Iffll The (l-R)V^ 
values are obtained using eq. 13.21 the experimental val- 
ues of v Q and 7oo0 for Cu +2 =-17.4 and Cu +1 =-5.2 in 
a crystalline environment. 

The calculated values produce the correct EFG prin- 
cipal axes. Furthermore, small values of the asymmetry 
parameters for Cu(2F/E) indicate a small deviation of 
these sites from tetragonal symmetry in accordance with 
our experimental results. The directions of the principal 
axes and the values of -q are not expected to be affected 
so much by the electronic structure of the subject ion 
since this effect is isotropic to first order. This explains 
the success of the point charge model in predicting the 
principal axes and t] values. A larger value of rf at is ob- 
tained for the Cu(lF) site. This value is not as large as 
the observed one since it does not contain the asymmetric 
ligand contribution. 

Since the Cu(2F/E) and Cu(lF) sites have +2 va- 
lence, one would expect |45j a valence contribution, (1- 
R-)V™ ( « -70 MHz, to the EFG from a hole in the 3d 9 
state and a zero contribution for Cu(lE) sites in the 3d 10 
state. Our values obtained for this contribution to the 
EFG for Cu(lE) site is small and supports this charge 
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configuration for this site. We find a difference in va- 
lence contribution of 2.19 MHz between the Cu(2F) and 
Cu(2E) sites. Comparing this value to the experimentally 
observed 3.15 MHz, we conclude that the valence (lo- 
cal moment) contribution to the EFG difference at these 
sites is dominant. This can be understood as resulting 
from a different charge transfer for these two sites; the 
full chains act as the charge reservoir whereas the empty 
chains do not. Charge transfered from the full chains re- 
sides at both the Cu(2F) and Cu(2E) sites but with a 
larger amount at Cu(2F), the site closer to the full chain. 
The inadequacy of the point charge model is apparent for 
the Cu(lF) site not only in its prediction of the asym- 
metry parameter but also in the valence contribution . 
In this case asymmetric ligand effects are mostly respon- 
sible. 

To investigate the effects of the domain walls on the 
EFG values, we have introduced different types of do- 
main walls and then calculated the EFG parameters for 
different sites. This calculation shows that the domain 
boundaries affect the Cu(lF/E) sites the most and the 
Cu(2F/E) sites the least. Our results also show that the 
bulk values of the lattice contribution to vq are essen- 
tially recovered just one lattice constant away from the 
domain walls. 



the stoichiometric YBC06.5 than previously available. 
Knight shift and spin-lattice relaxation measurements 
are in progress to investigate the pseudogap. In addition, 
as mentioned above, the long chains appear to support 
Friedel-like oscillations nucleated at the chain ends 32] . 

Another consequence of oxygen ordering in the chain 
layer is the effect it has on the Cu02 planes via charge 
transfer. This is manifested in the Cu(2E/F) splitting. 
Our analysis concludes that most of this splitting is due 
to differential charge transfer from the chains. 

Finally, previous reports of static moments in 
YBC06.5 are inconsistent with our results. 
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APPENDIX 

A straightforward calculation^^ to second order re- 
veals the transitions frequencies for spin 1=1 : 



IV. CONCLUSIONS 

Line assignments for YBC06.5 have been made. The 
lack of an observable NMR signal from the chain ends 
at the domain walls confirms the long chain lengths sug- 
gested by X-ray diffraction. We place a lower limit on 
the chain length of 50 lattice constants. Consistent with 
this is the narrow linewidths comparable to those of the 
best YBC07 crystals. This is a result of little oxygen 
disorder in these crystals. 

Thus these crystals are a much better realization of 



l/ = „(°) + „U + „W (A.l) 
^ (0) =j n H (l + K) (A.2) 
z/ 1 ) = — i>q(m — — ) \icos 2 9 — 1 — r/sin 2 8cos2(j)\ (A. 3) 



,( 2 ) 



57 2 
{102m(m - 1) - — } cos 2 9(1 + -rjcos2(j)) 



j n H (l + K) { 32 

9 2 
- {6m(m - 1) - -}(1 - -7jcos2(j)) 

n 2 r 21 
+ ^ |24m(m - 1) - 6 - {30m(m - 1) - —}cos 2 t 



(A.4) 



.51 



51. 



{— m(m - 1) - — }cos 2 2(/)( 



If 



where K is the Knight shift obtained by 

K = K x sin 2 9cos 2 4> + Ky sin 2 9sin 2 <j) + K z cos 2 9 (A.5) 

and 9 and (j) are the polar angles of the applied magnetic 
field with respect to the principal axis frame which is 



assumed to be the same for both the Knight shift tensor 
and EFG tensor. j n is the nuclear gyromagnetic ratio 
and H is the applied magnetic filed. ?y is the asymmetry 
parameter, defined as r\ — ss^Ym ; w here \V XX \ < \V yy \ < 
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\V ZZ \ are the principal components of the electric field 
gradient (EFG) tensor at the Cu site and vq = eQV zz /h, 
e the elementary electric charge, and Q is the electric 



quadrupole moment of the nucleus under study. The 
central and satellite transition resonance frequencies are 
obtained for ra=\ and m=|, —5, respectively. 
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